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Using the adiabatic connection, we formulate the free energy in terms of the correlation function of 
a fictitious system, h\{ r, r'), where A determines the interaction strength. To obtain h\(r, r') we use 
the Ornstein-Zernike equation, and the two equations constitute a general liquid-state framework 
for treating inhomogeneous fluids. As the two equations do not form a closed set, an approximate 
closure relation is required and it determines a type of an approximation. In the present work we 
investigate the random phase approximation (RPA) closure. We determine that this approximation 
is identical to the variational Gaussian approximation derived within the framework of the field- 
theory. We then apply our generalized RPA approximation to the Gaussian core model and Coulomb 
charges. 

PACS numbers: 


I. INTRODUCTION 

Pair interactions of hard-sphere fluids derive from 
the excluded volume effects: non-overlapping configura¬ 
tions recover an ideal-gas behavior, but the exclusion of 
overlapping configurations reduces available phase-space, 
leading at high density to phase transition. In this sense 
the hard-sphere fluids constitute a geometric problem. 
Within various successful (nonlocal) density functional 
theories (DFT), a free energy functional for hard-sphere 
fluids is built from a weighted rather than local density 
non-locality is attained by construction [l]. In early 
prescriptions, a weighted density corresponded to a con¬ 
voluted local density, where the single convoluting func¬ 
tion was the Mayer f-function. The resulting theories, 
their refinements and extensions came to be known as 
the weighted DFT theories. A crucial next development 
was to decompose a Mayer f-function into several weight 
functions, yielding multiple weighted densities and, by 
the same token, muiltiple building blocks from which an 
approximate F ex was to be constructed Referred 

to as the fundamental measure (FM) DFT, a nice feature 
of this approach is the capture of a correct dimensional 
crossover: each consecutive reduction of the system di¬ 
mensionality, 3D—>-2D—UD—>-0D, recovers either an ac¬ 
curate or exact F ex . 

The success of the hard-sphere DFT theories (and the 
lack of equivalent theories for arbitrary pair interactions), 
prompted attempts to implement the hard-sphere frame¬ 
work to other types of short-range interactions. It be¬ 
came something of a standard method to map parti¬ 
cles with arbitrary short-range interactions onto a hard- 
sphere fluid by ascribing to a pair potential of interest 
an effective diameter. Density profiles are then obtained 
from the DFT theories for hard-spheres. The Barker- 
Henderson effective diameter is one recipe among others 


for extracting an effective diameter [?}. 

A more sophisticated example is the ’’soft” funda¬ 
mental measure DFT developed for penetrable spheres 
(spheres may overlap but at an energy cost). Within this 
method F ex is constructed to satisfy a correct dimen¬ 
sional crossover [§l-IT2j. 

But for particles with arbitrary pair interactions, 
where excluded volume effects are not dominant, the 
mean-field approximation is still a preferred theoretical 
tool [IM3- An obvious example are charged particles 
with long-range interactions. Other examples are parti¬ 
cles with bound (non-divergent) interactions, known as 
penetrable particles. This class includes the Gaussian 
core model, or the already mentioned penetrable spheres. 
But the correlations neglected by the mean-field descrip¬ 
tion are not always trivial. This is particularly true of 
Coulomb systems. In such a case the ” beyond-mean- 
field” approach is desirable. Splitting the excess free 
energy into the mean-field and correlation contribution, 
fex[p] = \ f dr f dr' p(r)p(r')u(r, r') + F c [p), where p( r) 
is a number density, and u( r, r') is an arbitrary pair in¬ 
teraction, the ” beyond-mean-field” approach amounts to 
finding an appropriate functional F c [p\. 

The ” beyond-mean-held” approximations for Coulomb 
systems are dominantly formulated within the field- 
theoretical framework based on mathematical transfor¬ 
mation of a partition function, using a Gaussian integral 
identity [l8j . into a functional integral over an auxiliary 
fluctuating held [H-H- The saddle-point of the effec¬ 
tive Hamiltonian recovers the mean-held solution, while 
the harmonic fluctuations around the saddle-point ac¬ 
count for weak (Gaussian) correlations. If formulated 
variationally, the equations become self-consistent (non- 
perturbative) and generally are deemed superior to the 
perturbative formulation |l8j, HI, HU • 

A drawback of the held-theoretical formulation is the 
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loss of physical intuition after one moves from physical to 
auxiliary phase-space. In the present work we re-derive 
the variational Gaussian equations of the field-theoretical 
framework using only the liquid-state theory. The Gaus¬ 
sian approximation is equivalent to the well established 
random phase approximation (RPA) (27i|. a mathemati¬ 
cal signature of which is its b eing comprised of an infinite 
summation of ring diagrams [281 1 • Our formulation of the 
RPA is general and in principle applicable to any pair in¬ 
teractions. We apply our generalized RPA approximation 
to the Gaussian core model, considered to be a weakly 
correlated fluid [HI, a one-component plasma, and finally 
a symmetric electrolyte. 

In Sec. El we formulate the free energy within the 
liquid-state formalism using the adiabatic connection. 
By coupling it to the OZ equation, we set up a general 
theoretical framework for inhomogeneous fluids. In Sec. 
ED we consider the RPA closure and derive the appro¬ 
priate self-consistent equations. In Sec. IIVI we generalize 
the equations to multiple species. Finally, in Secs. El 
ED and ED we test the RPA approximation on concrete 
systems with wall geometry. 


II. ADIABATIC CONNECTION 
FORMULATION OF THE FREE ENERGY 


respectively. The free energy of a physical system can be 
expressed in terms of a fictitious system, as a thermody¬ 
namic integration, 

F = F ° + l! dX Trt ■ (5) 

where the reference free energy is 

F 0 [p\ = F id [p) + J dr p{r)U\—o{r), (6) 


and 

Fid\p] = k bT j dr p(r) [ log/?(r)A 3 - l] (7) 


is an ideal-gas free energy. The integrand in Eq. © can 
be written as 


dF x 

dX 


dr p(r) 


dU x (r) 

dX 


\ J dr J dr ' P(r)p(r')u{r,r') 

\ j dr j dr' p(r)p(r')h A (r, r')u(r,r') 


( 8 ) 


Given a general Hamiltonian for a system of interacting 
particles, 


and Eq. © becomes 

F[p\ = F id [p]+ J drp(r)U x = 0 (r) 


N - N 

H = '52U(r i ) + -J2u(. r i’ T j), (!) 

*=i 

where U(r) is an external potential, u( r,r') is a pair in¬ 
teraction, and N is the number of particles, our aim is 
to obtain a free energy expression in terms of physically 
meaningul quantities. To this end we use the adiabatic 
connection route dHH, wherein interactions are grad¬ 
ually switched on within a fictitious A-dependent system, 


+ J dr J dr' p(r)p(r')u(r,r') 

+ \ J dr f dr'p { r)p { r>Mr,r') J dXh\(r,r'). 

( 9 ) 


Then after a few cancelations we arrive at the final form, 


N . N 

H x = ^2 u x (n) + -^u( r i, r j), (2) 

*=i 

where the A-dependent external potential U x ( r) is in¬ 
troduced to keep a density fixed at its physical value 
for all A, and A = 1 recovers the physical potential, 
f/ A =i(r) = C/(r). 

The partition function and the free energy of a ficti¬ 
tious system are 

Z x = J dr 1 ... J dr N e~ 0Hx (3) 


F[p\ = F id [p]+ J drp(r)U(r) 

+ J dr J dr' p(r)p(r')u(r,r') 

+ i J dr J dr' p(r)p(r')u(r, r') J dXh x (r,r'). 

( 10 ) 

Note that the final expression does not depend on the 
fictitious potential U x . The only quantity that depends 
on A is a correlation function h x ( r,r / ). The last line of 
the expression represents the correlation free energy, 


and 


Fc[p] 



dr p( r) 


dr' p(r')h x (r, r')u(r, r') . 

(A) 


I3F X = log Z a , 


( 4 ) 
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Not surprisingly, F c depends on the correlation func¬ 
tion, h\( r, r'), that is obtained from the Ornstein-Zernike 
equation (OZ), 

h\{ r, r') = CA(r, r') + f dr" p(r")h\{r', r")cA(r, r"), 

. . ( 12 ) 

which is a well known exact relation within the liquid- 
state theory. Because the direct correlation function 
c\(r, r') is not known, an appropriate closure relation 
is still required. 


III. RANDOM PHASE APPROXIMATION 

We consider the simplest closure available, 

C x( r i r ') = -W r ,r'), (13) 

known as the random phase approximation (RPA). The 
closure modifies the exact Ornstein-Zernike relation in 
Eq. unj, 

M r > r ') = —/3Au(r,r') 

— /3A f dr" p(r")h x (r', r")u(r, r"). 

(14) 

Accordingly, we refer to it as the OZ-RPA equation. Ap¬ 
plication of the OZ-RPA modifies the correlation free en¬ 
ergy in Eq. CD, 

= <Mr) f\dX-^ J **)• 

(15) 

In the above equation zt(0) = u(r,r). Likewise, for a 
homogeneous fluids in a bulk we write hb( 0) = hb( r, r). 


A. Connection with the field-theoretical 
formulation 

The A-dependence in F c can be eliminated by expand¬ 
ing h x { r,r'), 

h x (r,r') = -/3Xu(r,r') 

+ (3 2 A 2 * * * J dr\ p(ri)u(r, ri)u(ri, r 7 ) 

- /3 3 A 3 Jdr iJ dr 2 p(ri)p(r 2 )u(r, ri)u(ri, r 2 )it(r 2 , r') 

+ ... (16) 

The expansion is generated iteratively by repeated inser¬ 
tion of the right hand side of Eq. 0 for every occur¬ 
rence of h x (r , r'). The notation is simplified by introduc¬ 
ing an operator 

, ( 17 ) 


and adopting a convention 

A n = JdriJdr 2 .. .jdr n - 1 A(r,r 1 )A(r 1 ,r 2 ) ■.. A(r„_i,r') 

(18) 

by means of which we get 


p(r)h x (r, r') = -\A + A 2 A 2 - A 3 A 3 + ... 
XA 

I+ XA 


(19) 


where / = A(r, r') is the identity matrix in the continuum 
limit. Integration over A now is done explicitly, 


" * 1 p(r)h x (r, r') . A 2 A 3 


dA ry ~ r AV 
A 


= -A+ —- —+■■■ = -1og[I+A], 

( 20 ) 


and F c [p] becomes 

F c [p] = -^Tr log[J + A] - ^ J dr p(r), (21) 

where the first term yields an infinite series of ring dia¬ 
grams, a characteristic feature of the RPA. 

The expression can further be rearranged by using the 
formal matrix identity, 


Tr log[7 + A}= log yj det [I + A] , (22) 


and the fact that a functional determinant is a solution 
of a Gaussian functional integral, 


1 


= / V(j>e 


~ \ I dr f dr ' 0(r)0(r')[i5(r”r')+A(r,r')] 


\J { det [/ + A] 

(23) 

where </>( r) is a fluctuating field and f V<j) is a functional 
integral. The partition function within the RPA can now 
be written as a Gaussian functional integral, 


Z rpa = e^“(°) e -^ mf 


x / Vtpe 2 


i / dr / dr' 0(r)0(r')[d(r-r')+A(r,r')] 


(24) 


where we used F = F m f + F c and Z = e~^ Fmt e~^ Fc . 
The functional integral formulation has been recovered 
without resorting to the Hubbard-Stratonovich transfor¬ 
mation, starting from the liquid-state formulation. 


B. density profile 

To obtain an equilibrium density we use the known 
thermodynamic condition, 


SF 

Sp(r) 


= /b 


A(r, r') = j3p(r)u(r, r') 


(25) 
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where f.i denotes the chemical potential. The functional 
derivative of F c with respect to p( r) incidentally elimi¬ 
nates all A-dependence, 


SF C 

Sp{r) 


k B T 6Tr log [/ + A] 1 

-3 — m — 2" (0) 


1 

2 L 


u(0) + k B Th( r, r) 


(26) 


and the functional derivative is written in terms of a cor¬ 
relation function of a physical system, h( r, r). The num¬ 
ber density that results is 

= /C , h g-/5C ; '(r) e -/3/dr'p(r')«(r,r') e i[/3M(0)+?i(r,r)]+/3/2 ex ^ 

(27) 

where we separated a chemical potential into ideal and 
excess parts, p = p t d + p e x, with the ideal contribution 
related to a bulk density, 


Here we assume that each particle feels different external 
potential, and pair interactions between different pairs 
are different. Of course, particles are not all different 
but are grouped into species. 

Assuming K different species, the free energy from adi¬ 
abatic connection is 

F [{Pk}} = F id [{p k }] +Y^ fdrU k (r)p k (r) 
k-r' 

+ \Y^j dr J dr ' Pk{r)pi{r')uu(r, r') 

+ dr J dr> Pk(r)pi{r')u k i(r,r') j^dXh^( r, r'), 


k,l 


where the ideal-gas contribution is 


(34) 


/ e ^id \ 

(28) 

K r 


V A 3 j- 

P F id[{Pk}} = J dr Pk ( r ) 

i= 1 J 

log Pfc( r )A 3 - 1 


Pb = 


The excess chemical potential within the present approx¬ 
imation is 


If the Ornstein-Zernike equation for multiple-species is 

K 


P ex = Pb [ dru(r) ~^[u(0) + k B Th b {0)], (29) h^(r,r') = c^(r,r')+5Z/ rfr Vn( r '0^( r V'04n( r V), 

J ^ n=l d 


where h b {r) is a correlation function in a bulk. For 
U (r) = 0, we accurately recover a bulk density, p( r) —> 
p b . More conveniently, a density can be written as 

= p foe —-St/( r ) e —/3 / dr' (p(r')-pb)u(r,r') e | [/»(r,r)-/i 6 (0)] _ 

(30) 


C. pressure 

Another quantity of interest is pressure that can be 
obtained from a type of thermodynamic integration in¬ 
volving a chemical potential, 


(36) 

where correlations between particles of a species k and 
l are mediated by all particles disregarding their type, 
then the RPA closure, = (3\uki(r, r'), yields 

hki( r,r') = -/3\u k i(r,r') 

~ I dr "Pn(r")hni('r\r")ukn(r , \r), 

n ^ 

(37) 

and the RPA correlation free energy is 

1 J'" f 

K pa = / dr ^( r ) 

/_1 j 


P — 

x ex — 


r pb , dp, 
d PP 


dp 


(31) 


k =1 
K 


- / drp k (r) 


k= 1 ‘ 


h kk ( r > r ) 

dX ~ 0A - ‘ 


(38) 


where P = k B T p b + P ex . The resulting expression shows 
A-dependence, 


j /* 7 rj~\ ” n 1 

P ex = t } pI J dru{r) -1 —p b h b { 0) - dX 


d b, a(0) 


(32) 


IV. MULTIPLE SPECIES 


We next generalize the RPA to multiple species. The 
fictitious Hamiltonian, equivalent to that in Eq. 0, is 


N 


A 


N 


H\ = U i( r i) + 9 5 Z U d( r i’ r j)- 


(33) 




(Compare with Eq. (11511 for a one component system). 
The lack of dependence on inter-species correlations, that 
is, hki{ r, r') for k ^ l, at first glance appears inaccurate. 
But as correlations between particles of the same species, 
hkk{ r, r'), are mediated by all the particles disregarding 
their type, the cross-correlations are always implicit in 
hkk (U r'). 


A. density 

An equilibrium density of a species k is obtained from 
the condition 

SF 

= Pk, (39) 


Spk( r) 
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and the correlational counterpart yields 

SF C u kk (0) k B T 

= - (40) 

If the excess chemical potential of a specie k is 

r ip i 

/4 X = E^ / dru k i{r)--\u k k(0) + k B Th b kk (0)\, (41) 

i J 

then a density is 

_ p\ e ~PU k {r) e -P Ef = i / dr' (p, (r')-pf ) Uw ( r ,r') 
x e 5 D l fcfc( r > r )—^fcfcCo)] _ M 2 ) 


where pb = TT 3 ^ 2 a 3 pb is the reduced density and V is the 
volume of a system. The correlation free energy density, 
f c = Fc/ V, becomes 

ep b ^ (—e ??) n_1 
/c 2 ^ n 5 / 2 

n=2 

= _aJ 1 + iMz£!?n (47) 

2 1 £ 1 ] J 

where Li m (;c) = is a polylogarithm. We may 

now obtain any quantity of interest. For example, the 
excess chemical potential is 


o £ L . Li 3/2 [~£V] 

^"■ = £ ’’-2\ 1 + - Tr, - 


(48) 


B. pressure 


Comparing with Eq. Cffl) we get another useful quantity, 


To obtain the pressure we use P ex = Y^ k =i PkPT ~ /ex, 
where the excess free energy density in a bulk is 


* 6 ( 0 ) 


F vH)>r 

£ n 3/2 


(49) 


/ex = ^J2p b kP b i I dru k i(r) 
k,i J 

- )i>£U(0) + t B r/Z^ 


fe=l 

The excess pressure then becomes 
i r 

Pex = ^PkPl / dlCU M( r ) 
1 J 


k, 

k B T 


k f z-i 

E^ Mo)- / dA 

k=i L • y ° 


4>) 


(43) 


(44) 


V. THE GAUSSIAN CORE MODEL (GCM) 

We apply the developed RPA approximation to the 
Gaussian core model (GCM), whose pair interactions 
have the Gaussian functional form, 

f)u{r) = ee~ r ! a . (45) 

a is the length scale that determines the interaction range 
and e determines the interaction strength. Because the 
potential is bound, the GCM particles are said to be 
penetrable. 

For a homogenous system the free energy in Eq. m 
can be calculated exactly, and each individual ring term 
becomes 


and the pressure is written as 

/3Pex £ f , Li 5 / 2 (-£r?) -U 3/2 (-eri) 

p b 2 { ep 

Our primary interest, however, lies in the RPA as a 
theory of inhomogeneous fluids. Considering a fluid con¬ 
fined by a hard wall at x = 0 to a half space x > 0 , we 
can use the contact value theorem to predict the density 
at a contact with a wall from a bulk pressure, 



p(o) = PP 

= Pb 


' f »7 Li 5 / 2 (—£ 77 ) - Li 3 / 2 (—er?) 

+ 2 + 2p 


(51) 


The first term is the ideal-gas contribution, the second is 
the mean-field contribution, and the last term accounts 
for the RPA correlations. 

In Fig. 0 we compare the contact density at a wall 
as a function of £ for different approximations. The 
RPA correlations lower the mean-field predictions, and 
the RPA corrections become more accurate at high den¬ 
sities, that is, for a larger number of overlaps. 

In Fig. © we plot the entire density profiles near a 
planar wall at x = 0. The mean-held is not expected to 
be accurate for e > 1 and the largest deviations from the 
exact results occur near the wall. The RPA profile not 
only improves the contact region but an entire profile, 
even for as large values of the interaction strength as 
£ = 7. 

To complete the analysis, we consider next a two com¬ 
ponent GCM system with interactions 


Tr A n 


V{£Vb) n 

(n7T<7 2 ) 3 / 2 


(46) 


Uij{x) 


2 / 2 

£e r 

2 / 2 

—£e r 


if i = j 
if i^j. 
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FIG. 1: The contact density, r/(0) = n 3 ^ 2 a 3 p(0), as a func¬ 
tion of an interaction strength e for a one component GCM. 
r/i, is the bulk reduced density. The dotted horizontal line 
corresponds to an ideal-gas prediction. 


The bulk density of both species is the same, Pb = pt = 
Pb , so that the mean-field contributions are canceled out 
and the density profile is determined strictly by correla¬ 
tions. In Fig. d3j) we compare a density profile of the 
RPA approximation with that from the simulation. The 
depletion of particles from the interface region is caused 
by unfavorable energy cost when a particle is removed 
from a bulk, which requires breaking of various ’’bonds” 
with its neighbors. 


VI. ONE-COMPONENT PLASMA 

We consider next Coulomb particles and transform the 
previously obtained expressions of the RPA approxima¬ 
tion to more familiar expressions in terms of an electro¬ 
static potential. The resulting expressions are the same 
as those obtained for a variational Gaussian approxima¬ 
tion within the field-theoretical framework [23|, H 3 and 
without resorting to the Hubbard-Stratonovich transfor¬ 
mation. 

Coulomb charges q interact via the following pair po- 



FIG. 2: Density profile for the GCM near a planar wall for 
different values of the interaction strength: e = 5, 6 , 7. A sim¬ 
ulation box dimensions are 20 a: 20 a: 20 a and it contains 1000 
particles. Without hard wall constraint this corresponds to 
a bulk reduced density r/b ~ 0.7. The numerical data points 
for the mean-field and the RPA correspond to the same con¬ 
ditions. 



FIG. 3: The total density profile near a planar wall for the 
two component GCM fluid, for e = 4. The system size and 
the number of particles is the same as in Fig. ©. except the 
total density is made up of two species which have identical 
density profiles. The density is determined by correlations 
alone as the mean-field effects cancel out. 

tential, 


where e is the background dielectric constant. A number 
density of Coulomb charges, using Eq. (13(J1) . is 

p{ r) = p 6 e-^ r M [M r , r )-M 0 )] ? ( 53 ) 

where the external potential, in electrostatic problems 
generated by permanent charges distributed over surfaces 
and accounted for by the boundary conditions, is omit¬ 
ted from the expression. Furthermore, we introduce an 
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electrostatic potential, ?/>(r), defined as 

qtj)( r) = J dr' p(r')u(r,r'). (54) 

To transform the OZ-RPA equation in Eq. (fl4l) into de¬ 
sired form, we apply the Laplacian operator to both sides 
of the equation, 


^72./ / \ fiq 2 
V h( r,r ) =- 


S (r - r') + p(r)h(r, r') 


where we used the identity 

V 2 n(r,r') = -(^-) d ( r-r'). 

We carry out the same operation on Eq. ([54]) . 

eVV(r) = -qp{ r), 


(55) 


(56) 


(57) 


and the result is the standard Poisson equation. 

Eq. (1531) . (l55l) . and (l57l) constitute the RPA approxi¬ 
mation for a density distribution and an electrostatic po¬ 
tential within the RPA level of approximation. Correla¬ 
tional contributions enter through the correlations in the 
number density h(r, r'). These can be related to correla¬ 
tions in electrostatic potential using a slightly rearranged 
OZ-RPA equation, 


h( r, r') = 


where —dU{z)/dz = —ecr c /2 is a constant force felt by 
particles on account of a uniform wall charge. The mean- 
field solution to this problem is 


Pmf{x) = — 


1 


1 2 


1 + (3qcr c x/2e 


(61) 


and it captures a weakly correlated limit. On the oppo¬ 
site end is the strong-coupling limit j31j, ill], 


p sc ( x ) = M e -/W/2^ (62) 

As correlations become significant, the density evolves 
from one functional form to another, p m f —> p sc . A per¬ 
turbative Gaussian approach for a counterion-only sys¬ 
tem yields a semi-analytic expression for a density cor¬ 
rection A p(z), p{z) = p m f(z) + Ap(z) [H, [32|, where 
f Q dz A p(z) = 0, to maintain neutrality, and Ap(0) = 0, 
not to violate the contact value theorem. The corrected 
density develops a ’’hump” at a short distance from a 
wall. It is difficult to justify or trace the hump to a 
physical cause as it is not confirmed by simulations; sim¬ 
ulations always yield a non-monotonic density profile. It 
should be concluded that the evolution of the ” hump” is 
an artifact of the approximation. It is hoped that the 
self-consistent equations of the RPA method eliminate 
the ’’hump” in iterative steps and in closer agreement 
with the true system. Our computations, however, indi¬ 
cate that self-consistency only slightly alters the results 

of tll6 pc'T’t.nrHa t.i'up a'n'nrnftrVi ( q pp ln'icr 0 TVip mnr lu- 


—/3 / dr" 


p(r")h(r', r") + <5(r 7 — r") 


4 ( r ! r”)- 


(58) 


The term in square brackets, p(r")h{r",r') 1 is identified 
as a correlation hole generated by a fixed particle at r' 
and the delta function denotes the density of a fixed par¬ 
ticle. The integral on the right hand side can be rein- 
terpretted as a perturbation of an electrostatic potential, 
T (r, r'), caused by a fixed particle at r' (the total electro¬ 
static potential is i/;(r)-t-\P(r, r')). The OZ-RPA equation 
simply becomes 


h(r,r') = — /3q^f(r, r'), (59) 

and the proportionality between the two fluctuating 
quantities is established. Note that this is not an exact 
equality but a result specific of the RPA approximation. 

For a concrete example we consider a counterion only 
system confined to a half-space x > 0. The counterion 
charge is q = e. The wall surface charge cr c at x = 0 
assures neutrality of the system. The dielectric constant 
is the same on both sides of the wall. As the bulk den¬ 
sity far away from the wall vanishes, the contact density 
is determined solely by the surface charge (not the pres¬ 
sure), 


f°°, , ,dU(z) Pal 

p(o) = ~J q dz P^— = ^ 


( 60 ) 




FIG. 4: Counterion density profiles for a counterion-only wall 
model. E denotes the coupling constant corresponding to the 
ratio of the Bjerrum and the Gouy-Chapman length, As and 
1/(27 tAbct c ), respectively. The larger the E, the more signifi¬ 
cant the correlations. 

sion is that the RPA approximation is not very accurate 
for the counterion-only system, and modifies the pertur¬ 
bative results only negligibly. 





















VII. ELECTROLYTE 


Continuing with electrostatics, we move toward elec¬ 
trolytes. We consider a symmetric electrolyte, q : q, with 
a bulk concentration of both species p The two types 
of pair interactions are 


J rt(r, r') if k = 

( —u( r, r') if k^l 


where u{ r, r') is the Coulomb potential given in Eq. (f52)l . 

For a two species system, there are four different cor¬ 
relation functions, hki( r, r'). Within the RPA they can 
be expressed in terms of a single function. Accordingly, 
we have 


as was previously demonstrated for a one-component 
plasma in Eq. ( 1551 ) . 

Note that the excess chemical potential does not in¬ 
clude the mean-field contributions and depends exclu¬ 
sively on correlations, 


P 


± 

ex 


1 

2 


lim 

r —>0 


kBTh b (r) + u(r) 


(69) 


The mean-field contributions cancel out by virtue of 
charge neutrality. The same is true of pressure which 
reads 


li- 


ksTpb 

2 


lim 

r —>-0 




h b( r ) 

A 


(70) 


, / / h(r, r') if k = l 

h ki( T > r ) - | -h(r,r') if k^l 

and h( r, r') is obtained from the OZ-RPA relation, 

h{ r,r') = 


-P / dr" 


p(r")h(r', r") + S(r' — r") 


u(r,r"), 


(63) 


where p( r) = p+(r) + p-(r) is the total density. Within 
the RPA, the number density of each species is 

P± (r) = p b e^^ei [m^-MO)] ) ( 64 ) 


where the correlation function is obtained from a trans¬ 
formed Eq. (1551) (by applying the Laplacian operator to 
both sides of the OZ-RPA equation), 


eV 2 /i(r, r') = (3q 2 


p(r)h(r, r') + S(r 



(65) 


Together with the Poisson equation, 

eV 2 ip( r) = -p c (r), ( 66 ) 

where p c ( r) = qp+(r) — qp~{ r) is a charge density, we 
have a complete approximation for a density and electro¬ 
static potential of a symmetric electrolyte. 

As in the case of a one-component plasma, we may link 
the correlations in a number density to the correlations in 
an electrostatic potential. We identify the term in brack¬ 
ets in Eq. (l63l) as a charge correlation hole generated by 
fixing either a positive or a negative charge g at r'. For 
a fixed positive charge we have 

Phoie(r,r') = p+(r)/i ++ (r,r') - p_(r)/i+_(r,r') 

= p+(r)h(r,r') + p_(r)h(r,r') 

= p(r)h(r,r'). (67) 


Consequently, the fluctuations in the number density are 
proportional to the fluctuations in electrostatic potential, 


( 68 ) 


and where p b = p b + + p b _ is the bulk total density. 

For bulk electrolytes Eq. (l65l) recovers the Debye- 
Hiickel theory for a point-charge, 


d 2 hj(r) 

dr 2 



(71) 


where k\ = \J\/3q 2 pb/e is the screening parameter. The 
Debye-Huckel solution is 

A f3q 2 e~ Kxr 


and the excess pressure becomes 


(72) 


p — 

1 ex 


q 2 Pb 


lim 


87 re 0 

k q 2 Pb 


-I dX ■ 


3 8 « (73) 

where k = K\—i . The total pressure may be written as 


^< 74 » 

We note that the RPA correlations reduce the ideal gas 
pressure. From the contact value theorem we may infer 
that at neutral interfaces the contact density will be lower 
than that in a bulk, indicating a depletion zone that is 
not caused by dielectric discontinuity but a more efficient 
bonding arrangement within a bulk. 

Once again resort to a simple wall geometry for test¬ 
ing purposes. The wall is uncharged and its only func¬ 
tion is to create an interface between an electrolyte and 
an empty space. Dielectric constant is uniform across 
the interface and everywhere else. The results in Fig. 
(Taj) indicate the depletion zone near an interface, neither 
caused by dielectric discontinuity nor finite ion size but 
correlations. The repulsion of ions from an interface in 
turn increases the surface tension [35[. As for the two 
component GCM in Fig. m this is caused by a more 
efficient salvation of ions in a bulk, where each ion is ac¬ 
companied by an opposite charge-cloud. Near the inter¬ 
face, the charge-cloud is deformed by a nearby interface 
rising the cost of an energy. 


h(r,r') = -/3g^(r,r'), 






















9 



FIG. 5: Density profiles for a 1:1 electrolyte near a wall 
at z = 0. The Bjerrum and the screening lengths are 
Ab = 0.72 nm and k - 1 = 0.33 nm, respectively. The ions 
are depleted from a wall region to minimize the number of 
particles at an interface. 

VIII. CONCLUSION 

The free energy formulation using the adiabatic con¬ 
nection offers a simple way for incorporating the OZ 
equation into a theoretical framework for liquids and 
may be considered as an extension of the integral equa¬ 
tion theories to inhomogenous fluids. Finally, the choice 
of closure determines an approximation. In the present 
work we explore the RPA closure, CA(r, r') = — X/3u(r, r'). 
We demonstrate that the resulting general RPA approx¬ 
imation is equivalent to the field-theoretical variational 
Gaussian approximation, derived completely within the 
liquid-state theory. 

We test the developed RPA framework for different in¬ 
homogeneous fluids for a simple wall geometry. For the 
Gaussian core model the RPA density profiles show de¬ 
cisive improvement over the mean-field and are in good 
agreement with a simulation, even for large interaction 
strengths. For the counterion-only system the RPA is 


less accurate. Like the perturbative Gaussian approxi¬ 
mation [22| [33] the density profile develops an unphysical 
bump near the wall. In comparison, the simulated pro¬ 
files are always monotonically decreasing. Furthermore, 
the self-consistency of our approach appears to modify 
rather negligibly the profile obtained from the pertur¬ 
bative scheme [22|, [33[ . Consequently, we conclude that 
the RPA (or the variational Gaussian approximation) is 
not an accurate theoretical tool for strongly correlated 
Coulomb fluids. Finally, we test apply the RPA to a 
symmetric electrolyte near a neutral interface, without 
dielectric discontinuity. We observe the depletion of den¬ 
sity near an interface generated exclusively by correla¬ 
tions, since the mean-field contributions in this system 
are cancelled out. 

As final remarks, we restate that we did not see that 
the self-consistent scheme of the RPA (or the variational 
Gaussian) produces significant modifications in compari¬ 
son with the results obtained perturbatively. If there are 
situations where self-consistency is crucial, we cannot be 
sure, but for systems and parameters considered in this 
work we did not come across such conditions. Finally, as 
the future project, we think it worthwhile to explore the 
adiabatic connection framework presented in this work 
but for more accurate closures. Self-consistency in these 
more advanced closures may turn out to be more signifi¬ 
cant. 
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